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The combination of current large scale structure and cosmic microwave background (CMB) 
anisotropics data can place strong constraints on the sum of the neutrino masses. Here we show that 
future cosmic shear experiments, in combination with CMB constraints, can provide the statistical 
accuracy required to answer questions about differences in the mass of individual neutrino species. 
Allowing for the possibility that masses are non-degenerate we combine Fisher matrix forecasts for 
a weak lensing survey like Euclid with those for the forthcoming Planck experiment. Under the 
assumption that neutrino mass splitting is described by a normal hierarchy we find that the combi- 
nation Planck and Euclid will possibly reach enough sensitivity to put a constraint on the mass of 
a single species. Using a Bayesian evidence calculation we find that such future experiments could 
provide strong evidence for either a normal or an inverted neutrino hierachy. Finally we show that 
if a particular neutrino hierachy is assumed then this could bias cosmological parameter constraints, 
for example the dark energy equation of state parameter, by > lcr, and the sum of masses by 2.3a. 



I. INTRODUCTION 

Accurately determining the absolute value of the 
neutrino mass is one of the main goals of particle 
physics. However since the effect of individual neu- 
trinos is small it is cosmological observations, that 
observe the cumulative effect of neutrinos on large 
scales, that present the most powerful way to bound 
the absolute neutrino mass scale. Albeit indirect and 
model dependent, cosmological constraints are cur- 
rently stronger than those coming from beta-decay 
experiments (for recent reviews see Jl[ and Q). For 
example, the Mainz 0] and Troitsk [J Tritium decay 
experiments give upper limits on the single electron 
neutrino mass of m < 2.05eV and m < 2.3eV re- 
spectively, at 95%c.L. Cosmological data coming from 
Cosmic Microwave Background (CMB) measurements 
from WMAP [f| experiment combined with baryonic 
oscillation data Q and various Supernovae observa- 
tions (0, H) lower this limit to to < 0.22eV 0. On 
the other hand observations of flavour oscillations in 
atmospheric and solar neutrinos provide evidence not 
only for a non-zero neutrino mass, but also for a dif- 
ference between masses, measuring squared mass dif- 
ferences of fTo| : 
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where ranges indicated are 3a confidence level and 
toi, TO2 and to 3 are the three mass eigenstates. The 
ambiguity in the sign of Ato 2 ^ leads to an uncertainty 
about the neutrino mass scheme, allowing for two pos- 
sible hierarchies: the normal hierarchy, given by the 
scheme to 3 3> TO2 > toi, or the inverted hierarchy 
to 2 > mi 3> to 3 . Note that given equation (fl|) an in- 
verted hierarchy scenario would be automatically ex- 
cluded by measuring a total mass X) 771 " ~ O.leV. 



It is commonly perceived that cosmology is able to 
constrain the total neutrino mass ^ rriy while mass 
differences between eigenstates can be neglected. This 
is an excellent approximation as shown in 11| at least 
for the CMB anisotropics power spectrum, since the 
effect of neutrino mass on the CMB is related to the 
physical density of massive neutrinos uj v , i.e. to their 
total mass ^2 m y . Individual neutrino masses do have 
an effect on the matter power spectrum, due to the 
different transition redshifts from relativistic to non- 
relativistic behaviour. This effect is still much smaller 
than that due to the total mass and can be safely ne- 
glected in analysing current cosmological data. Nev- 
ertheless in the near future various experiments will 
reach a much higher accuracy in reconstructing the 
matter power spectrum. It is therefore timely to con- 
sider the possibility that these surveys will be sensitive 
to single neutrino masses. In recent papers (see for ex- 
ample [l2j and [HI]) a forecast has been made consid- 
ering future observations of CMB anisotropics, CMB 
lensing and galaxy distribution finding that this kind 
of data doesn't seem able to reach enough accuracy to 
discriminate between the two hierarchies. In [l4| has 
been shown that future accurate measurements of the 
redshifted 21 cm signal from the epoch of reionization 
can in principle measure individual masses but will 
be very difficult to achieve the precision required to 
distinguish normal and inverted hierarchies. 

In [l5j an explicit and more general parameteriza- 
tion of neutrino mass splitting was introduced, rep- 
resenting various possible hierarchies; with the only 
simplifying approximation that two neutrinos are of 
the same mass m\ = to,2- This approximation is 
well justified by equation ((T|). In [15| it was shown 
that even a cosmic variance limited CMB experiment 
would not be able to detect a difference in individual 
neutrino masses. In addition CMB lensing, even if 
limited by cosmic variance only, is strongly inhibited 
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in measuring the mass hierarchy by degeneracies with 
other cosmological parameters. 

In this article we use a Fisher matrix formalism ap- 
plied to the same parameterization of [l5[ to assess the 
ability of future cosmic shear measurements, like those 
achievable with Euclid/DUNE [lj] experiment, com- 
bined with Planck [ItJ CMB data to place constraints 
on single neutrino masses. It is almost a decade that 
cosmic shear has been recognized as one of the most 
powerful tools to constrain the total neutrino mass 
(see for example [lU, [IH and[2(|) and hence an inves- 
tigation into how far these experiments can proceed 
in the exploration of the neutrino properties is well 
justified. 

The article is organized as follows: in section UT1 we 
review in more detail the effects of the total neutrino 
mass and of individual neutrino masses on cosmolog- 
ical observables. We also describe the parameteriza- 
tion of [l5[ which will be used throughout the rest of 
the article. In section HTT1 we describe the implemen- 
taton of our cosmic shear and CMB Fisher matrices. 
Section IIVI shows results from our forecasts as func- 
tion of various assumptions about the neutrino mass 
hierarchy, and for different parameter sets. We also 
discuss results in the light of our parameterization of 
the neutrino mass hierarchy. Our forecasts show that 
accurate measurements of the matter power spectrum 
from Euclid, combined with CMB data from Planck, 
can reach the accuracy required to constrain single 
neutrino masses. In section [V] we conduct a more ac- 
curate analysis of our results through a Bayesian ev- 
idence investigation. We also show that assuming a 
wrong hierarchy can lead to a bias in the recovered 
value of other cosmological parameters, in particular 
for the dark energy equation of state, generally com- 
parable with the lcr sensitivity. The largest bias is 
found in the total neutrino mass, due to the degener- 
acy involving ^2 m v and the hierarchy parameter, at 
2.3cr. Finally, in section fVTl we summarize our confu- 
sions. 
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II. EFFECTS OF NEUTRINO MASS 
SPLITTING 



The main effect of a non-zero neutrino mass on 
cosmology is through the collisionlcss fluid behaviour 
that causes neutrinos to free-stream over a typical 
length which is known as free- streaming length Xfs- 
The consequence of this free propagation is a can- 
cellation of neutrino density fluctuations on scales 
smaller than Xps an d a slow-down in the growth 
of perturbation on these scales. The matter power 
spectrum results are then damped for wave vectors 
k > kps — 2tt/Xfs- The free streaming wave vector 
of a single species Icfs depends on the mass of that 
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where h = Hq/ (100km s 1 Mpc 1 ). On these small 
scales the matter power spectrum is suppressed with 
respect to the power spectrum of a cosmology with 
masslcss neutrinos by an amount that depends mainly 
on the fraction /„ of matter density in the form of 
massive neutrinos (/„ = Q u /Q m ). In the case of de- 
generate masses f2„ can be expressed as a function of 
the total neutrino mass: 
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However, as pointed out in [22j, even in the case of 
non-degenerate masses this relation remains a good 
approximation. 

Massive neutrinos become non-relativistic at a red- 
shift given by z nr ~ 2 • 10 3 m y /eV, so that neutrinos 
with masses up to ~ 0.5eV are still relativistic at time 
of recombination. As a result the effect of neutrino 
free streaming on the CMB power spectrum is negli- 
gible for small neutrino masses. In this case the main 
effect of neutrino masses on the CMB is indirect, re- 
lated to the delay of matter-radiation equality. This 
causes a small shift in the peaks of the power spec- 
trum and a slight increase of their heights due to a 
longer duration of the Sachs- Wolfe effect. 

From this discussion is clear that neutrino mass 
affects growth of structure in two ways: the matter 
power spectrum is suppressed by an amount that de- 
pends mainly on the total mass ^ m v but also, even 
if for a minor amount, on single masses, because the 
time of transition to non-relativistic regime depends 
on single masses. The typical wave vector over which 
this suppression can be observed also depends on the 
mass of single neutrino species as shown by equation 
([2|). Therefore a reconstruction of the matter power 
spectrum can in principle give information both on 

m v an d on single masses (through the suppression 
of the power spectrum and &fs), even if this second 
effect is generally much smaller than that due to the 
total mass. In the following sections we investigate the 
ability of future cosmological experiments to reach the 
sensitivity required to detect differences in neutrino 
masses. Following fl5| we parameterize mass splitting 
by introducing the parameter a defined as the fraction 
of the total mass in the third neutrino mass eigenstate: 
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The other two eigenstates are assumed to share the 
same mass mi = m.2- This approximation is sup- 
ported by the observed differences in the squares of 
neutrino masses as measured by oscillations of atmo- 
spheric neutrinos, |Amf 3 |/Amf 2 ~ 0.5-10~ 2 [Io| . The 
advantage of this parameterization is that it allows us 
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to represent in a simple way various mass hierarchies 
such as total degeneracy (a = 1/3) and normal or 
inverted hierarchy given respectively by a ~ 1 and 
a < 1. 



The Fisher matrix for weak lensing is given by (e.g. 
0) 

WL v (2l+l)Al ! 5P fcm ! 



III. 



FORECAST FOR WEAK LENSING 
TOMOGRAPHY 



Weak lensing (see 2JJ for a recent review or 
http://www.gravitationallensing.net) is a par- 
ticularly powerful probe of cosmology since it simul- 
taneously measures the growth of structure through 
the matter power spectrum, and the geometry of the 
Universe through the lensing effect. Since weak lens- 
ing probes the dark matter power spectrum directly 
it is not limited by assumptions about galaxy bias 
(how galaxies are clustered with respect to the dark 
matter). Future weak lensing surveys will measure 
photometric redshifts of billions of galaxies allow- 
ing the possibility of 3D weak lensing analysis (e.g. 
[32 . [33l . [Hi], Hfl) or a tomographic reconstruction of 
growth of structures as a function of time through a 
binning of the redshift distribution of galaxies D(z), 
with a considerable gain of cosmological information 
(e.g. on neutrin os jl3| : dark energy [35|; the growth 
of structure (37l . l38l| and map the dark matter distri- 
bution as a function of redshift [39|). 

In this section we present forecasts for the upcom- 
ing weak lensing survey Euclid [l6|, combined with 
constraints expected from the CMB Planck experi- 
ment [l7j using a Fisher matrix formalism. We first 
explored a simple 8-parameter model with fiducial val- 
ues fl b fi 2 = 0.022, Q c h 2 = 0.111, h = 0.7, r = 0.084, 
n s = 0.95, As = 2.48 • 10" 9 , £m„ = 0.055eV and 
the hierarchy parameter a — 0.95. Note that we are 
assuming a fiducial normal hierarchy scheme for neu- 
trino masses. We then repeated the analysis for a 
larger set of parameters, including running of spec- 
tral index (with target value dn s /dhik = 0) and the 
dark energy equation of state w (assuming cosmolog- 
ical constant w = — 1 as the fiducial model), to study 
how degeneracies with various parameters affect con- 
straints on neutrino mass hierarchy. 

To calculate the power spectra we used the numeri- 
cal code CAMB [23| that can calculate CMB and mat- 
ter power spectra, also for non-degenerate neutrino 
masses. The usual definition of the Fisher matrix [30| 
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where L is the likelihood function for a set of param- 
eters pi. When the derivatives of ([!]) are evaluated at 
the fiducial model the Fisher matrix gives an estimate 
of the best statistical error achievable on the param- 
eters (via the Cramer-Rao inequality) for the method 
and survey design considcrd 



where Pij{£) is the convergence weak lensing power 
spectrum that depends on the non-linear matter 
power spectrum at redshift z, P n i(k,z), obtained by 
correcting the linear matter power spectrum P{k, z) 
using the option halof it of CAMB A£ is the 

step used for I and 
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In the last expression 7j„ t is the rms intrinsic galaxy 
ellipticity (and we assume { / jf nt ) 1 ^ 2 = 0.22) and 
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the number of galaxies per steradian belonging to j th 
bin 
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where d is the number of galaxies per square ar- 
cminute and hj is the fraction of sources belong- 
ing to the j th bin. For the Euclid experiment we 
take d = 35 and f s k y = 0.5. The galaxy redshift 
distribution is assumed to have the form D(z) oc 
z 2 exp[— (z/zo) 1 - 5 } with z = 0.9. For this experiment 
photometric redshift uncertainties are assumed to be 
a z = 0.03(1 + z). We treat this uncertainty follow- 
ing the approach of [26| where a galaxy with red- 
shift z could be wrongly observed at a redshift z p h- 
Letting p(z\z p h) be the probability that this happens 
the distribution of galaxies in the i th bin is modi- 
fied to take into account this 'leakage' between bins 
Di(z) = J z h . dz p hD{z)p(z\zph) ■ We choose a simple 
Gaussian form for p{z\z p y l ) 



p(z\z ph ) = (2na z ) 
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The Fisher matrix for weak lensing is then added 
to that of CMB to obtain constraints from the com- 
bination Planck+Euclid 
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For a CMB experiment the Fisher matrix is given by 
El: 
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where the couples PP' and QQ' mean in our case 
TT, TE or EE (temperature and E-mode polarisa- 
tion). Cove is the power spectrum covariance matrix 
at the £ th multipole and f max is the maximum multi- 
pole available given the angular resolution of the con- 
sidered experiment, for Planck we use £ max = 2000. 
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The other specifications we used are listed in Table 
HI The total Fisher matrix F tot is then inverted to 
obtain uncerainties on cosmological parameters of our 
model. 



PLANCK Channel/GHz FWHM AT/T AP/T 
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TABLE I: specifications for Planck experiment used in the 
Fisher matrix calculation. AT/T and AP/T are sensi- 
tivities (fiK/K) for temperature and polarization respec- 
tively. 



IV. RESULTS 

In this section we show results of our forecasts on 
neutrino mass parameters ^ m v and a. 



A. Parameter Constraints 

The constraints of Table [TT| and Fig. [1] show that 
for ^2 niv = 0.055eV, if the neutrino mass hierarchy is 
described by a normal hierarchy (a = 0.95) then the 
combination of weak lensing data from an experiment 
like Euclid and accurate measurements of the CMB 
power spectrum, achievable with Planck, can become 
sensitive to the mass of single species. The 1<t un- 
certainties on the sum of masses is cr^p m„ = 0.037eV, 
which is in agreement with other forecasts (see for ex- 
ample [2(|) and confirms the ability of these surveys 
to detect neutrino mass. In particular this work has 
a resonance with Debono et al. (2009) (in prepara- 
tion) in which the effect of parameter sets on neutrino 
mass constraints are investigated. We find that for 
the parameters that are common between the two ar- 
ticles there is agreement between the predicted errors, 
despite the slightly different parameter sets, power 
spestrum estimation approaches (CAMB vs Eisen- 
stein & Hu code [2l[) and assumptions. 

As one can see for a typical normal hierarchy sce- 
nario with a = 0.95 the combination Euclid+Planck 
can reach a ~ 20% sensitivity on a with an error 
a a = 0.19. For this normal hierarchy scenario we 
have repeated the Fisher matrix calculation using a 
larger parameter space, including running dn s /dlnk 
of the spectral index and the dark energy equation of 
state parameter w, to check the weakening in the con- 
straints on a induced by degeneracies among these 
parameters. One may expect that running and w 
would have a large degeneracy with neutrino mass 
since they both can add an effective damping on small 
scale. However, as shown in in Table [TTI and in Fig. [TJ 



constraints on a are not seriously weakened - for this 
second 10-parameter space a a = 0.22. 

experiment cry- mv a a 

Planck 0.49 eV > 1 

Planck+Euclid 0.037 eV 0.22 

TABLE II: contraints from Planck+Euclid Fisher matrix 
on neutrino mass parameters for the 10-parameter space 
described in the text. 
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FIG. 1: 68% and 95% probability contours (two- 
parameter) in the plane m v f° r Planck+Euclid from 
our Fisher matrix calculation for the 10-parameter space 
described in the text. 



B. Parameter isat ion Investigation 

For the smaller 8-parametcr set we repeated the 
Fisher matrix calculation for different fiducial values 
for a. In Fig. [5] we show the relative uncertain- 
ties in a from the combination Euclid+Planck as a 
function of the target model. The figure is a com- 
bination of a general trend, that causes the uncer- 
tainty on a to decrease for a increasing, and the way 
in which constraints from CMB combine with those 
from weak lensing. The loss of sensitivity around 
a ~ 0.88 is in fact due to a strong rotation of the 
degeneracy in the plane m„ in the Euclid Fisher 
matrix, as shown in Fig. [31 from an anti-correlation 
to a positive correlation. For the Planck Fisher ma- 
trix instead the degeneracy between these two param- 
eters is less dependent on the fiducial value of a. As 
one can see for a ~ 0.88 the constraints from Euclid 
are almost independent of a, and the combination 
of the two experiments lose the ability to break the 
degeneracy between these parameters. This is con- 
firmed also by the fact that the only other parameter 
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that shows a significant peak of uncertainty in corre- 
spondence of a ~ 0.88 is the total mass ^2 m v itself 
that icrease from cv TOj/ — 0.025eV for a = 0.86 to 
cr E „ iiy ~ 0.034eV for a = 0.88. 

We emphasise that this is a peculiar effect of this 
parameterisation, we have tested the derivative nu- 
merically for convergence. Therefore we advocate this 
parameterisation but with a strong warning that re- 
sults arc highly dependant on the fiducial value of a, 
particularly around a m 0.88. 

We note however that the general decrease in <r a /a 
is quite intuitive since for 7773 ~ ^2 m v a variation in 
a means a variation only in the mass that can have 
an effect on the growth of structure. In this limit the 
others two eigenstates have a very low mass and are 
rclativistic up to a very low redshift. For example for 
our target model {^2,m v ~ 0.055cV and a = 0.95) 
the two eigenstates have mi = m 2 ~ 1.4 • 10~ 3 eV, 
this results in the neutrinos being relativisitc up to a 
redshift z ~ 2. Hence for a — ► 1 the sensitivity to 
this parameter increases. Conversely when a is sig- 
nificantly different from 1, for example in the case of 
total degeneracy, when a = 1/3, a variation in a and 
hence in 777,3 implies an opposite variation in mi and 
7772 which now have a non- negligible mass. These vari- 
ations partially compensate reducing the sensitivity of 
cosmology to a. 

This reasoning is confirmed by the comparison be- 
tween derivatives of the matter power spectrum P(k) 
(that enters in the calculation of the weak lensing con- 
vergence power spectrum) with respect to a calculated 
for different values of a. In Fig. [4] are shown deriva- 
tives dP/da for the case a = 0.95 and a — 1/3. The 
latter is about two orders of magnitude smaller than 
the first. As one can see the two derivatives have oppo- 
site signs due to the different change in kps induced 
by a change in a in the two target models. When 
a — ► 1 an increase in a causes an increase in fc^s ac- 
cording to equation ([2]) while, as we have said above, 
mi and m 2 are not important. Instead, for a — 1/3 an 
increase in this parameter causes a decrease of mi and 
mi (which now cannot be neglected) and the overall 
effect is a decrease in kps- 



V. BAYESIAN ANALYSIS 

To assess better the power of these cosmological 
probes to detect a neutrino mass difference we present 
a Baycsian evidence forecast. Calculation of Baycsian 
evidence allows one to make a comparison of different 
models, as opposed to parameter estimation within a 
model. According to Bayes theorem the (posterior) 
probability of a set of par ameters describing a model 
\i given the data d is [3l| 
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FIG. 2: relative error on a as a function of the target 
model for the combination Planck+Euclid. The fiducial 
value for the total mass is ^ m v — 0.055eV. 
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FIG. 3: 68% 2-parameter probability contours in the plane 
X] m v -a for Planck (green, light ellipse), Euclid (blue, dark 
ellipse) and the combination (red, central ellipse) from 
our Fisher matrix calculations for the 8-parameter space 
described in the text. The plots are for fiducial values 
a = 0.86 (left), a = 0.88 (middle) and a = 0.9 (right). 
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FIG. 4: derivatives of matter power spectrum with respect 
to q for two target models. The derivative for a = 1/3 is 
multiplied for a factor 10. See text for more details. 
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where C(9) = p(d\9,/i) is the likelihood function, 
p(9\fi) is the prior probability on the parameters of 
the model [i and p(d\p) is the Bayesian evidence; 
that has the role of a normalization constant, being 
p{d\n)=jd6p{d\6,n)p{6\n). 

The following is a summary of the technique de- 
scribe in (29|. Given a cosmological model M' de- 
scribed by a number of parameters n' , a common 
problem is to verify whether data require the inclusion 
of some new parameters in the model so as creating 
a new more complicated model M with a number of 
parameter n > n'. In this case one has to take the 
ratio of the posterior probabilities of the two mod- 
els p(M'\d)/p(M\d). This ratio can be obtained from 
Bayes thorem: 

p(M'\d) _ p(M') 

p{M\d) p{M) 1 ' 

where p(M') and p(M) are the prior probabilities of 
the two models and B is the Bayes factor given by 
the ratio between Bayesian evidences: 

Jd9'p(d\9',M')p(9'\M>) 
Jd9p(d\9,M)p(9\M) [ ' 

where 9' and 9 are the set of parameters of M' and 
M . In the case we are considering the two models 
are nested, in the sense that they share the same n' 
parameters. For nested models and for Gaussian like- 
lihoods (fTTjl approximates to [29[ 

B = ( 27 r)-P/ 2 ^p2J cxp (~69 a F a0 S9 ) f[ A9 n , +q , 
Vdet F' V 2 / ^ 

. . (12) 

where the index q runs over the p = n — n' additional 
parameters of M with respect to M' and A9 are the 
prior ranges on the parameters. Note that the Fisher 
matrix F is n x n while F' is n' x n' . 

The p parameters are assumed to be fixed at a cer- 
tain fiducial values in the model M' are shifted by an 
amount St/} with respect to their fiducial value in M. 
For a, /3 = 1, n' then this shift 89 a is given by 

56 a = -{F'- 1 ) a pG M 5^ 1 y=l,...,p, (13) 

where for a, (3 = l,...,p we have S9 a = Sip a . The 
quantity G that appears in (fT5| is a block n' x p of 
the full Fisher matrix F. 

The calculation of Bayes factor through Fisher ma- 
trices helps to clarify whether future experimental 
data will be sensitive to a wrong assumption about 
some parameters (for example fixing p parameters to 
wrong fiducial values). Bayesian analysis is known 
to be conservative in the sense that models with a 
smaller number of parameters are favoured until data 
strongly require the introduction of new parameters 
f3lj ] . Hence, for the case of nested models only very 
sensitive experiments will have the power to discern 



that a certain parameter is kept fixed to an incorrect 
value. 

In what follows we apply equation (|12[) to our 
neutrino Fisher matrices so that we can understand 
if future weak lensing and CMB data will achieve 
enough sensitivity to require a parameterization of 
the mass hierarchy. All results showed below are 
for the combination of the Euclid and Planck experi- 
ments and hence the Fisher matrix used in the cal- 
culation of Bayes factor is the total Fisher matrix 
F l ° = F^ L +F^ MB . We assume that the true model 
M is represented by the 10-parameter model described 
in the previous section with a normal hierarchy scheme 
for neutrinos (a = 0.95) and m v = 0.055eV. We 
next consider a simpler 9-parameter model M' in 
which the parameter a is fixed to a certain value, 
shifted of an amount Sa with respect to a = 0.95 
and calculate the Bayes factor for these two compet- 
ing models. 

Results of the evidence calculation are shown in 
Fig. [5] where we plot |lnS| as a function of \8a\. 
The horizontal lines corresponds to the values of the 
Jeffreys scale (36j : | \nB\ < 1 means inconclusive ev- 
idence, 1 < |ln_B| < 2.5 is a substantial evidence, 
2.5 < I hiB\ < 5 is a strong evidence and | ln_B| > 5 
is considered decisive evidence. The area leftward of 
the cusp corresponds to a Bayes factor B > 1 which 
in our case means evidence for the simpler model M' . 
For \8a\ > 0.24, B becomes smaller than 1 indicat- 
ing that data would require the introduction of a in 
the analysis, and will be sensitive to differences be- 
tween neutrino masses. In particular for da ~ —0.62, 
that represents the common (and under our hypoth- 
esis wrong) assumption of total degeneracy between 
neutrino masses, the data would give a strong evi- 
dence for model M requiring a parameterization of 
neutrino mass differences. 

If data will be able to give evidence for a hierarchy 
of neutrino masses then it is proper to verify what the 
effect of assuming degeneracy of masses, or a wrong hi- 
erarchy, on other cosmological parameters will be. As 
shown in p9j fixing one parameter to a wrong value 
causes a shift in the best fit value of other parameters 
according to (fl"3)) . We calculate the bias in cosmolog- 
ical parameters due to a wrong assumption for a and 
assuming normal hierarchy (a = 0.95) as true model. 
The results are shown in Table IIHI and show that as- 
suming a degenerate hierarchy (Sa = —0.62) or an 
inverted hierarchy (Sa = —0.9) would cause a shift in 
other cosmological parameters comparable with the 
la statistical error. Only the shift in ^2 mv is sig- 
nificantly greater than the la error due to the high 
degeneracy between a and ^ m v . 

These results confirm those of Fig. [5] the data will 
be accurate enough to require a parameterization of 
the mass splitting, this causes a non-negligible bias 
in other parameters in the case of a wrong assump- 
tion. However the bias on cosmological parameters 
is generally smaller than the ltr uncertainties. No- 
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FIG. 5: absolute value of In B as a function of \5a\. The 
lines indicates the limits of the Jeffreys scale. On the right 
of the cusp is B < 1, meaning evidence favours a more 
general parameterization of neutrino mass hierarchy. 
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FIG. 6: Jeffreys scale contours of In B as a function of 
8a and 8^2 m v . The inner part of the plot (from the 
innermost \nB = 1 contour) corresponds to values B > 1 
and hence evidence for the simpler model (see text). The 
cross indicates the assumption of degenerate masses fixing 
the total mass to the correct value; the square indicates a 
typical inverted hierarchy scenario. 



table exceptions include the dark energy equation of 
state parameter. Assuming total degeneracy of neu- 
trino masses, the shift on w is smaller than la but 
becomes slightly greater assuming an inverted hierar- 
chy. Our Fisher matrix analysis indicates that this 
is due essentially to the shift in ^2 n^v and hence to 
the degeneracy w- J] m„ because there is no significant 
correlation between w and a. 

We also show the results of an evidence calcula- 
tion for an even simpler model M", in which both 
the total mass ^ m v and a are fixed, over the model 
M. Note that now we are comparing an 8-parameter 
model with a 10-parameter model. The results are 







Inverted 


Degenerate 


Parameter 


a 


8a = -0.90 


8a = -0.62 


w 


0.041 


-0.047 


-0.033 


Q. b h 2 


10~ 4 


-0.4 • 10~ 4 


-2.7- 10~ 5 


Q c h 2 


0.00065 


0.0013 


0.00090 


h 


0.013 


0.0049 


0.0036 


T 


0.0028 


-0.0012 


-0.00082 


n a 


0.0022 


-0.0036 


-0.0024 


A, 


1.44 • 10 -11 


5.75 • 10" 12 


3.94 ■ 10~ 12 


£m„(eV) 


0.037 


0.086 


0.060 


dn s /dink 


0.0031 


-0.0019 


-0.0012 



TABLE III: la errors on cosmological parameters and the 
bias (8) due to a wrong assumption in the neutrino mass 
hierarchy. 



shown in Fig. [6] where we have plotted the contour 
values of | ]nB\ as a function of 5a and 5^2 m„. The 
inner contours from 1 to 5 correspond to values B > 1, 
and so to an evidence which favours the simpler model, 
while the outer contours are relative to values B < 1 
and hence to evidence for the more complicated model 
M . The star in the plot refers to the wrong assump- 
tion of total degeneracy (but fixing the total mass to 
the true value <5^m„ = 0). As one can see, under 
these assumptions the data would favour the simpler 
model giving a substantial evidence for M" . 

The evidence for the simpler model is due to the 
smaller number of parameters of M" with respect to 
M (the Occams razor term, see [2!|) and also to the 
negative degeneracy between a and £ m v (given by 
the off-diagonal term of the inverted full Fisher matrix 

(^E^aM^W^Em.Em, - ~0.44). 

Because of this degeneracy there is a region of confu- 
sion in the plane of Fig. [5] in which a wrong assump- 
tion in a is compensated by a wrong, and opposite 
in sign, assumption for ^ m v , leading to evidence in 
favour of the simpler model A/". This explains why 
in Fig. O for a fixed 5a the the Bayes factor initially 
increases becoming greater than 1 for 5 ^2 m u increas- 
ing. Of course when 5 J2 m v becomes large enough 
the Bayes factor decreases because the data start to 
favour the true model M. 

If one assumes an inverted hierarchy (represented 
by the square in the plot), for example fixing the to- 
tal mass to the minimum value allowed for this hier- 
archy m v — 0.1 leV) and a to a low value (we take 
a = 0.05), one obtains a strong evidence for the more 
complicated model M. 

These evidence calculations indicate that a future 
weak lensing survey could become sensitive to the hi- 
erarchy of neutrino masses requiring a suitable param- 
eterization of mass splitting. We note however that 
even if a model assuming total degeneracy could be 
still favoured. This indicates that Planck and Euclid 
data could give strong evidence for the existence of a 
hierarchy among neutrino masses. 
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VI. CONCLUSIONS 

In this paper we have investigated the ability of 
future cosmic shear measurements, like those achiev- 
able with the proposed Euclid mission, to constrain 
differences in the mass of individual neutrino species. 
Using an explicit parameterization of neutrino mass 
splitting and a Fisher matrix formalism we have 
found that the combination Euclid+Planck will reach 
enough sensitivity to put constraints on the fraction 
of mass in the third neutrino mass cigenstatc (a) 
under the assumption that the neutrino mass scheme 
is described by a normal hierarchy, with fiducial value 
a = 0.95: a a — 0.22. We have also investigated the 
parameterization dependence of our results, repeating 
our forecasts calculation for different fiducial values 
of a and finding that constraints on a are generally 
decreasing for increasing a. We have found a loss of 
sensitivity around a ~ 0.88 due to a strong rotation 
of the degeneracy in the plane ^ m v -a in the Euclid 
Fisher matrix, suggesting a considerable dependence 
of the results on the fiducial value around this point. 
We have then studied more deeply the power of de- 
tecting neutrino mass splitting considering Bayesian 
evidence. We have found that future data from 
these experiments will provide strong evidence for 
a neutrino mass splitting against a model assuming 
degeneracy of masses. As a consequence we have 



also shown that assuming a wrong hierarchy can 
bias constraints on other cosmological parameters, 
in particular with those parameters involved in 
degeneracies with Y) m v and a. For the dark energy 
equation of state we found a bias greater than the 
la statistical error {cr w = 0.041 from Euclid+Planck) 
assuming an inverted hierarchy S w = —0.047 and 
comparable with la assuming degenerate masses 
5 W = -0.033. 

In conclusion we emphasize that, even if these 
constraints are strongly dependent on the parame- 
terization used, the possibility of having a splitting 
in neutrino masses cannot be neglected in analysing 
future data from Euclid- like experiments. A wrong 
assumption about neutrino mass hierarchies can in- 
deed cause non-negligible bias on other cosmological 
parameters and in particular on the dark energy 
equation of state. 
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